Bipolaron-SO(5) Non-Fermi Liquid in a Two-channel Anderson Model with 

Phonon-assisted Hybridizations 
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We analyze non-Fermi liquid (NFL) properties along a line of critical points in a two-channel 
Anderson model with phonon-assisted hybridizations. We succeed in identifying hidden nonmagnetic 
SO (5) degrees of freedom for the valence-fluctuation regime, and we analyze the model on the basis 
of boundary conformal field theory. We find that the NFL spectra along the critical line, which is 
the same as those in the two-channel Kondo model, can be alternatively derived by a fusion in the 
nonmagnetic SO(5) sector. The leading irrelevant operators near the NFL fixed points vary as a 
function of Coulomb repulsion U; operators in the spin sector dominate for large U, while those in the 
SO(5) sector dominate for small U, and we confirm this variation in our numerical renormalization 
group calculations. As a result, the thermodynamic singularity for small U differs from that of 
the conventional two-channel Kondo problem. In particular, the impurity contribution to specific 
heat is proportional to temperature and bipolaron fluctuations, which are coupled electron-phonon 
fluctuations, diverge logarithmically at low temperatures for small U. 

PACS numbers: 75.20. Hr, 74.25. Kc 



I. INTRODUCTION 

Strongly interacting electron-phonon systems have at- 
tracted much attention in condensed matter physics. Vi- 
brating ion oscillations in metal interact with conduc- 
tion electrons, leading to interesting low-energy phenom- 
ena such as superconductivity and various density-wave 
states. For about three decades, Kondo effects^ due to 
local ion oscillations have been studied intensively by var- 
ious authors^ - — Realization of heavy-fermion behaviors 
in filled-skutterudite SmOs4Sbi2 under high magnetic 
fields^ suggests that this compound is not a conventional 
heavy-fermion system due to magnetic Kondo effects. 1 
One of the distinct properties in the filled-skutterudite 
structure is that the Sm atom is located inside a large 
"cage" consisting of Sb atoms. Since the size of the cage 
is much larger than that of the Sm atom, the spacial de- 
pendence of the potential energy for the Sm atomic os- 
cillations is very shallow and the oscillations become an- 
harmonic. Materials that have similar cage-like structure 
such as clathrate compounds^ and /3-pyrochlore oxides^ 
have also been studied recently, partially due to the po- 
tential application as thermoelectric materials and obser- 
vation of strong coupling superconductivity mediated by 
strong anharmonic oscillations. 

As a pioneering work of Kondo physics in electron- 
phonon systems, Yu and Anderson investigated a sys- 
tem with one local phonon interacting with spinless 
two-channel conduction electrons^ Vladar and Zawad- 
owski investigated so called two-level systems where 
"two-level" represents two quasi-degenerate ion states in 
a double-well potential, coupled with two-channel con- 
duction electrons, and they proposed that it is possi- 
ble to realize the two-channel Kondo phenomenal in 
this system. There are many theories that investigate 
such Kondo physics due to ionic oscillations j 7 ' 8 ' 11 Re- 



cently the full phase diagram of a two-channel An- 
derson model with phonon-assisted hybridizations was 
clarified^— by using Wilson's numerical renormaliza- 
tion group2£ (NRG) method to analyze the Kondo effects 
in molecular systems^ and also in cage compounds such 
as fillcd-skuttcruditesi22r— A line of two-channel Kondo 
like fixed points was found from the weakly correlated 
regime to the Kondo regime in the phase diagram. The 
effects of anharmonicity in the local phonon potential 
were also investigated^ 

The purpose of this study is to clarify why the two 
completely different regimes, i.e., the Kondo and weak 
coupling regimes in the model investigated in Refs. 19- 
21, can be connected smoothly along the critical line, 
where the fixed point spectra and the quantum numbers 
characterizing each of the states are invariant*^— In 
the Kondo regime, it is natural to expect that a mag- 
netic two-channel Kondo model (2CKM) 18 describes low- 
energy properties of this system. For the weak coupling 
regime, however, it is unclear what is going on. One 
naive expectation is that some nonmagnetic degrees of 
freedom play important roles to realize the identical fixed 
point spectra and residual entropy^ ln-\/2 with those in 
the magnetic 2CKM. In this respect, one can imagine 
that nonmagnetic two-channel Kondo effects occur, as 
had been proposed in two-level systems^— ^ where an 
ion tunnels between two local potential minima, inter- 
acting with two-channel conduction electrons. However, 
this expectation is not supported by the asymptotically 
exact NRG resultsJ^2i 

In this paper, we will demonstrate that the non-Fermi 
liquid (NFL) of the 2CKM can be alternatively inter- 
preted as a nonmagnetic SO(5) NFL. This point of view 
can resolve the above questions and we can further pre- 
dict various crossover behaviors along the critical line. In 
Sec. 2, we will introduce the SO(5) degrees of freedom 
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and develop a critical theory along the line of the fixed 
points by using the non-Abelian bosonization method 
and boundary conformal field theory. Section 3 will be 
devoted to confirming the analytic results obtained in 
Sec. 2 by using Wilson's NRG method. Finally, in Sec. 
4, we will summarize the present results. 

II. BOUNDARY CONFORMAL FIELD 
THEORETICAL APPROACH 

A. Two-channel Anderson Model with 
Phonon-assisted Hybridizations 

1. Model 

In this paper, we investigate a two-channel Ander- 
son model with one-component directional local phonons 
that assist hybridization processes between a localized 
electron at the origin and conduction electrons^ This 
model has been investigated in the context of the Kondo 
effects in molecular systems while taking into account its 
vibrations.— Later, the same model was reanalyzed to 
investigate the Kondo physics in cage compounds that 
include magnetic ions in their cage structure! 20 ' 21 ! 24 

The Hamiltonian is 

H = ^ e q (n qsa + n qpa ) + U n fa - l) 2 

qa a 

+ Yl^( v o^ + Vixp ga ) + h.c.} + n^b, (i) 

go- 
where s qa .(p qrT ) is the conduction electron creation opera- 
tor with the radial wavenumber q, the spin a and s-wave 
(p-wave) symmetry around the origin. e q represents the 
energy dispersion of the conduction electrons, and we set 
the total volume to unity, /T is the creation operator of 
a localized orbital with the spin a and n qsa = s qcr s qrT , 
n qpcr = PqaPqa, and rif a = f£f a . U, V , and V\ repre- 
sent the Coulomb repulsion, the hybridization, and the 
phonon-assisted hybridization, respectively. The model 
([I]) is essentially the same as the model used in Ref. 19, 
in which the parameters satisfy Vr = Vl and A = 0. 
The dimensionless displacement of p-wave local phonons 
is indicated by x = b + b* with b^ being the phonon cre- 
ation operator and is the energy of the phonon. For 
simplicity, we restrict ourselves to the particle-hole sym- 
metric case, since the particle-hole asymmetry does not 
alter our main conclusion. 



2. Symmetry 

Before going into detailed analysis, we list here the 
symmetries of the present system. The Hamiltonian ([1]) 
has three symmetries. One is the spin SU(2) symmetry 
and another is the charge SU(2) symmetry. The last is Zi 
symmetry which is related to inversion symmetry at the 



origin. For each of the symmetries, there are conserved 
quantities. 

For the spin SU(2) symmetry, the total spin and its 
z-component are conserved. The total spin is given by 

Stot = S + ^S s (.x / ) + ^S p (x / ), (2) 
i i 

where S is the spin of the /-electron and S s (a:;)(Sp(x/)) 
is the spin of the s(p)-electron at site X\ in a one- 
dimensional "radial" lattice. Since S to t satisfies SU(2) 
commutation relations, the eigenvalue of S 2 ot = j(j + 1), 
with j being half integers or integers. 

With regard to the charge SU(2), it is well known that 
the total axial charge and its z-component are conserved. 
The axial charge for the /-electron I = (I x ,I y ,I z ) is given 
as 

I « = KS n /'- 1 ) 1 (3) 

a 

1+ = ll=l x +lly = -f\fl (4) 

The axial charges for the s and p-electrons are defined in 
the same way as S s ^(xi), and we denote them as l s (%i) 
and I p (xi), respectively^ Then, the total axial charge is 
given as 

Itot =l + Y, I ^ x ') + ^2M x i)- ( 5 ) 

The axial charge operators satisfy SU(2) commutation 
relations and thus, the eigenvalue of I 2 ot is i(i + 1) with 
i being half integers or integers. 

Finally, for the Zi symmetry, the total parity P is 
conserved: 

P = mod(^n gpCT + 6 t 6,2V (6) 

The eigenvalue of P is or 1. The values and 1 cor- 
respond to even and odd parity, respectively. Note that 
under the inversion operation, s qa — > s qa , p qa — > —p qa , 
and b — > —b. Here, since q is the radial wavenumber, q 
does not change. 

3. Background 

In this subsection, we explain what remains to be clar- 
ified in this model and what we have already understood 
in the previous worksJ^— 

As we have mentioned in Sec. 1, the global phase di- 
agram of the model ((T|) is known in the Vi-U plane. We 
schematically draw the phase diagram in Fig. [TJ There 
are two phases in the V\-U plane with fixed f2 and Vq£& 
Each phase is characterized by the ground state par- 
ity. In the phase for small Vi, it is even parity, while 
it is odd parity for large V±. For large U, there appears 
spin-1/2 local magnetic moment of the /-electron in both 



phases. The spin is eventually screened via Kondo cou- 
plings by that of the s- or p-electrons, as determined by 
the strength of Vb and V\. For large U, the phonon state 
can be regarded as an even-parity state denoted by |e) 
in Fig. [U which is continuously connected to the phonon 
vacuum state. Between the two phases, there is a line of 
NFL fixed points characterized by the spectra equal to 
those in the magnetic 2CKM. 

As U decreases, the local magnetic moment of the /- 
electron disappears and the Kondo-singlet state in the 
spin sector gradually crossovers to the different configura- 
tions in each of the phases. For small V\ and the small U 
regime, the ground state is essentially a non-interacting 
state, which was called "renormalized Fermi liquid" in 
Ref. 20. For large V\ and small U, the ground state con- 
sists of, in addition to the component dominant for large 
U, the odd-parity phonon state |o) coupled with even- 
parity states formed by /-and p-electrons. Note that 
near the critical line, there are components of spin-singlet 
states between / and p(s) electrons with \o) for the left 
(right) side of the line (not depicted in Fig. [l| and their 
magnitudes are comparable with those of the spin-singlet 
states with |e). 

A natural question is that, why the line of the NFL 
fixed point is continuous, even when the dominant com- 
ponents of the f -electron in the ground states crossover 
from magnetic to nonmagnetic ones. It should be also 
clarified that, why even when the local magnetic moment 
is absent for small U , the spectra of the NFL is the same 
as those in the magnetic 2CKM. 

In Sec. Ill Fit we will construct a critical theory that 
can describe the line of the NFL fixed points. Although 
it is beyond the scope of this paper to investigate physics 
very far from the critical line, it is possible to analyze 
the stability of the NFL fixed points and predict various 
critical behaviors in our theory. 



B. Non-Abelian Bosonization 




FIG. 1: (Color online) Schematic ground-state Vi-U phase 
diagram of model (fl) . 20 ' 21 There are two phases. One is the 
phase characterized by the parity-even ground state for small 
V\ , while the other is the odd-parity state for large Vi . Inside 
the each phase, the electron and also the phonon configura- 
tions crossover as U varies. |e)(|o)) represents the even (odd) 
parity phonon state. Between the two phases, there is a line 
of NFL fixed points for all the values of U > where the 
fixed point spectra of the NRG^ 2 ^ are the same as those in 
magnetic 2CKM. In the figure, we depict characteristic con- 
figurations with the spin singlet in each region in the two 
phases. f-s(p) indicates the spin-singlet state between / and 
s(p) electrons and ff(pp) represents the double occupied state 
for the f(p) electron. 



apply non-Abelian bosonization to this model, the sim- 
plest way of so called conformal embedding is to bosonize, 
global charge U(l), spin SU(2) and flavor SU(2) degrees 
of freedomi 30 i 31 We set the one-dimensional system size 
as [—1,1] and the bosonized Hamiltonian leads 



In this section, we map the original Hamiltonian to 
one in an effective one-dimensional continuous "radial" 
space with only left-moving conduction electron compo- 
nents, and we apply non-Abelian bosonization^ In this 
approach, the free electron part of the Hamiltonian (fTJ) 
is written as 



H = — - I dx 

zir 



sl(x)d x s r7 (x) + pl(x)d x p a {x) , (7) 



where x is the position in the one-dimensional space and 
vp is the Fermi velocity. 



1. Conformal embedding: U (1)®SU (2)i®S\J (2)i 

The Hamiltonian ([7]) consists of two flavors of conduc- 
tion electrons s and p with the spin j — 1/2. When we 



Hn = 



7TVF 



, ■ F n • F_ n 



(8) 



where J n , J„ = (J*, jv, J*), and F„ = (F£, F%, F*) rep- 
resent the charge, spin, and flavor (left moving) current 
operators in the Fourier space labeled by integers n, re- 
spectively, and : A : is the normal ordering of operator A. 
J n satisfies the U(l) boson commutation relation, and J„ 
and F„ satisfy the SU(2)fe Kac-Moody algebra with the 
level k = 2. This conformal embedding is suitable when 
the interactions consist of, for example, exchange interac- 
tions in the spin sector as in the case of the conventional 
two-channel Kondo problem . 30 ! 31 However, it is inconve- 
nient for us, since there is no flavor SU(2) symmetry in 
the Hamiltonian (fTJ). 
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2. Conformal embedding: SU (2)2® SO (5)2 

As is well known, the symmetry in the 2CKM is higher 
than U(l)(g>SU(2)<g>SU(2), and it is SU(2)®SO(5),22 In 
the NRG studies of the Hamiltonian Q, SU(2)®SO(5) 
symmetry is also realized along the line of the NFL fixed 
pointSi 20 ! 21 In this subsection, we clarify what SO (5) de- 
grees of freedom are in the model ([T]) . 

First, we introduce the Nambu representation: 

f(x) = [s\(x),- P \(x), Pi (x), Sl (x)}, (9) 
*(x) = 'Mx),-^),^),^)]. (10) 

We find that the SO(5) "density" is given by using ^(x) 
as 



4 4 

L ah (x) = ^2^ a (x)(C ab ) a ^ ? (x), 
0=1 0=1 



(11) 



where C ab = T ab /2 with 1 < a < b < 5 are SO(5) genera- 
tors and defined as 4 by 4 matrices as shown in Appendix 
They define SO(5) rotations and satisfy the SO(5) 
commutation relation (|A3|) . There are ten generators in 
the SO (5) group, which are the adjoint representation, 
denoted by 10. For later purposes, we define the 10- 
component vector L(x) as 

L(x) = (L 12 ,L 13 , i 14 , L 15 , L 23 , L 24 , L 25 , L 34 , L 35 , L 45 U2) 

where x dependence is omitted on the right-hand side 
of Eq. (fT2"j). We can bosonize L(x) and their Fourier 
components L„ satisfy the following SO(5)/c Kac-Moody 
algebra with the level k = 2^ 

i L n b ' L m} = l f ah ' Cd ' ef L E J +m + —5 a c5bd5 n +rn,a-{^'i) 

Here, f ab cd > e f is the SO (5) structure constant and is 
given by Eqs. (|A"2l and ([AS]) . 

The Hamiltonian ([7|) is bosonized by using the confor- 
mal embedding SU(2) 2 ®SO(5) 2 , leading to 



#0 - 12 [4 : J ™ ' J 



(14) 



Indeed, the central charge for the SO(5) sector is cgo(5) = 
5/2 and that for the spin sector is c s = 3/2, leading to 
c s + c,so{5) — 4 as it should be. This form (fl4|) is used 
in a spin-3/2 dipole-octupole Kondo problem.— There, 
\£(x) corresponds to the spin-3/2 fermion operator and 
the SO(5) generators correspond to linear combinations 
of the spin-3/2 dipole and octupole operators, and the 
spin current J is replaced by the SU(2) axial charge cur- 
rent. Since, in the spin-3/2 model, quadrupole operators 
are classified as the SO (5) vector, i.e., the 5 representa- 
tion, we define the corresponding degrees of freedom in 
our model and they are given by 



4 4 



where r a 's with 1 < a < 5 are 4 by 4 matrices defined in 
Appendix IA 21 In the following, we use a five-component 
vector n(x) = (n 1 (x),--- , n 5 (x)). Note that the spin 
eigenvalue of n(x) is j = 1 and the z-component j z — 0, 
while j ' = for L(x). j z = ±1 components of the 5 
representation n± (x) can be constructed by applying the 
spin operators J ± (x) = J x (x)±iJ v (x) to n(x). In total, 
28 operators, J, L, n, and n±, form a complete set of 
the conduction electron "density" operators in the sense 
of the Nambu representation. 

It is also important to check the "independence" of 
the two sectors and we find, by direct calculations, 
[L ab (x), J^(x')] = 0, i.e., they are "independent." The 
actual form of L(x) and n(x) is given in TablcUl It is clear 
that L(x) consists of charge, flavor, and spin-singlet pair- 
ing operators, while n(x) consists of spin-triplet pairing 
and a complex object of spin-flavor operators, see Ap- 
pendix [Bj In terms of the charge current J n and the 
flavor current F„ in Eq. §8§ are related to — i ? \ 5 , and 
{~Ll\Ll\Lf l respectively. 

By using this conformal embedding, the free electron 
energy spectra E are calculated via eigenvalues of the 
Casimir operators in the two sectors and are given by-2^ 



En — 



7TVF 



i(j' + l) , ^50(5) 



(16) 



(15) 



a=l 0=1 



Here, m is a non-negative integer, m > and Cso(5) 
is the eigenvalue of the Casimir operator in the SO (5) 
sector. The eigenvalue depends on the irreducible rep- 
resentations and is given as Cso(5) — for the identity 
1, 5/2 for the spinor 4, and 4 for the vector 5. Eigen- 
states with m — means that the eigenstate is a primary 
state, while those with m > 1 indicates that the states 
include particle-hole excitations and are classified as de- 
scendant states in the conformal tower characterized by 
a set of primary states^ The SU(2)2 Kac-Moody alge- 
bra restricts possible values of j for the primary states in 
the spin sector. The spin j of the primary states should 
be < j < k/2 = lM The SO(5) 2 Kac-Moody alge- 
bra also restricts the number of primary states in the 
SO(5) sector, and, as shown in Appendix O there are 
three primary states in the SO (5) sector: the identity 1, 
the spinor 4, and the vector 5. Using Eq. (HHJ), we can 
reproduce the free-electron spectra as shown in Table |TT] 
(a). 



3. Local operators classified in a hidden SO(5) group 

In this subsection, we will introduce local "flavor" de- 
grees of freedom. Using them combined with I, we will 
show that we can construct local SO(5) degrees of free- 
dom in terms of local operators: /, f\ b, and b^ for small 
U. 

In Ref. 21, it is shown that the energy spectra of 
small-cluster problems capture the essential aspect of the 
critical line obtained in the NRG calculations. Namely, 
there is a level crossing of the ground states of the 
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TABLE I: List of operators in the SO(5) sector. Each of 
irreducible representation is labeled by (jj z , dimension of the 
SO(5) group). Position index x is omitted for L(x) and n(x) 
and the operator forms shown are those multiplied by a factor 
2: 2L(x), 2n(x), 2Lj and 2nj. Each operator is labeled by the 
parity P, the axial charge eigenvalue i and the z-component 
i z , and the flavor F (see, analysis in Sec. Ill B 4|) . ii z — l±i 
means that the operator is in linear combinations of i z = 1 
and — 1. 
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Cm) = (0, 1/2) sector. Here, j(i) represents the eigen- 
value of total spin (axial charge). Let us briefly explain 
their results in the following. 

First, for large Vi, it is natural to consider a two-site 
problem, where / and p(0) electrons and b are taken into 
account. For [7 = 0, the Lang-Firsov transformation 36 
provides the exact solution of this problem. The ground 
state for U = is doubly degenerate. The two are two- 
electron states with the spin being singlet. This degen- 
eracy is distinguished by the parity, i.e., one is an even- 
parity state \4> e ) and the other is odd: \<j> ). This degen- 
eracy, however, is lifted by U, and for small but finite U, 
the ground state is the even-parity state with a small ex- 
citation gap to the odd-parity state. Importantly, in the 
low-energy states for U — 0, the phonon b appears only 



cosh 
sinh 



2A(6-& t )]|0) ph , 
2X(b - &t)l |0) ph . 



(17) 
(18) 



Here, |0) p h represents the vacuum of the phonon and A = 
—Vi/Cl. The other phonon states are in higher energy 
above f2, and thus do not play important roles. 

The three-site problem, where /, p(0), and s(0) elec- 
trons and the phonons b are taken into account, exhibits 
a qualitatively correct phase diagram, if we see the sec- 
tor of two-electron or four-electron states, i.e., i = 1/2 
states. For small V\, the ground state of the two-electron 
sector is an even-parity spin-singlet state, \<p' e ) ~ \<fie), 
while for large V\ the ground state changes to an odd- 
parity spin-singlet state \<p' ) ~ \4>o)- The same is true in 
the four-electron sector, and we denote them as \ip' e ) and 
\ip' D ). Their spin eigenvalues are also j — 0. Thus, there 
is a line where a level crossing occurs. Along the level- 
crossing line, the energies of these four states — \4>' e )i W a )i 
\t/j' e ), and \ip' ) — coincide and they form an SO(5) spinor, 
4 representation, in the SO(5) group. With regard to the 
= (0,1/2) sector, there are quasi-degenerate states 
with even and odd parity that originate in the ground 
states for the two-site problem with one electrons in the 
s-orbital. One is an SO(5)-singlet, while the other is one 
of the states in the 5 representation. 

Our assumption on the phonon degrees of freedom for 
small U is that the only two states, such as |e) and |o), 
which are solely constructed by the phonon part, are im- 
portant in the low-energy properties along the critical 
line. Indeed, this is valid in the limit of U = 0, since the 
first excited state lies at fi, and thus it does not play any 
role in the low-energy physics.^ Then, we can construct 
"flavor" operators characterized by the Pauli matrices 
in these two bases and we denote them as t x ,t v , and t z , 
identifying the even-parity state as Iff) and the odd-parity 
state as | JJ-) . The level crossing in the three-site problem 
is described by the presence of a term oc (U—U c )t z while 
fixing other parameters, where U c is the level-crossing 
value of U . Note that, for U > 0, there is always a fi- 
nite gap between the even and odd ground states in the 
(0,1/2) sector, since the gap arises mainly from the fact 
that the two are classified in a different irreducible rep- 
resentation and thus have generally different energies. 

Another local degree of freedom is the / operator part. 
We consider quadratic operators in terms of f a and /J. 
Possible linear combinations are the spin S and the axial 
charge I; see Eqs. © and (g]). Since the SO (5) 

sector is nonmagnetic, we consider the /-electron part of 
the axial charge, I, in detail. 

I satisfies the SU(2) commutation relations, and thus, 
the eigenvalue of I 2 is + with i = or 1/2. When I a 
acts on an i = subspace /||0), where |0) is the vacuum 
of the /-electrons, l a f^\0) — 0. Because of this property, 
the operation of I a is automatically projected on an i = 
1/2 subspace, |0) and /1/f |0). This is very important in 
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the derivation of local SO(5) degrees of freedom below, 
since the algebra of I a = 2I a is very similar to that of 
the Pauli matrices. Indeed, J a 's satisfy 







(19) 




- VP8 a p, 


(20) 


V 


= it = ly = Iz = 2n ft n fl ~ Iz, 


(21) 


Lv 


— VI — I 


(22) 



where V is the projection operator onto the i — 1/2 sub- 
space and is alternatively given by V = 4I 2 /3. 

Now, let us first introduce an SO (5) vector, i.e., the 5 
representation n/: 



n/ = ]:{IxTy,VT x ,I z Ty,VT z 



-h T v) 



(23) 



Note that n/ is a spin-singlet operator unlike n(x) of 
the conduction electron with the spin being triplet, and 
nj includes I a in all the components, which means n/ 
acts on the % = 1/2 subspace. At this stage, it is not 
proved that n/ transforms as an SO(5) vector, since we 
have not defined any local SO (5) generator. Thus, let us 
define a trial local SO(5) generators, 10, via Eq. (|A5|) . 
2Lf = [2n?,2n|]/(2i), and we obtain 



(24) 



The trial SO(5) generators Lf"s, indeed, satisfy 
the SO(5) commutation relations ([A3]> : [Lf,Lj d ] = 
-i^Lf-SacLf-SbdLf + SadLf). Thus, Lf denotes 
the SO(5) generators that we seek. We can also confirm 
that the commutation relation between 2L/ and 2n/ is 
similar to Eq. ([Al5]) . [2nJ,2iJ c ] = -i(sgn(c - a)5 a6 nj + 
sgn(& — a)5 ac rtj) i.e., n/ is the 5 representation of the 
local SO(5) group defined by L/. 
Finally, we note that the four states 



-f}flm,f}flm) 



(25) 



transform as an SO (5) spinor, 4 representation. This can 
be checked by applying ladder operators in the SO(5) 
group listed in Appendix [Cl 



4- Exchange interactions 

So far, we have analyzed how SO (5) degrees of free- 
dom appear in the free part of the Hamiltonian ([7]) and 
what local SO(5) operations for small U are. In this sub- 
section, we phenomenologically discuss possible and/or 
impossible interactions under the SU(2)<£>SO(5) symme- 
try, which are closely related to two fusions introduced in 
Sec. Ill B 51 We show two possible exchange interactions 
but they are not the "microscopic" effective Hamiltonian 
of ((T|). They should be regarded as one of the effective 
interactions for the coarse-grained system. 



The simplest invariant form under SU(2)®SO(5) sym- 
metry is exchange interactions between the impurity and 
the conduction electron component at the origin x = 0. 
Spin-spin exchange interactions are evidently possible 
under SU(2)® SO (5) symmetry: 



H S = JS- 3(0), 



(26) 



with J being the coupling constant. This term preserves 
the symmetry of the original Hamiltonian (JTJ) , and can 
be regarded as the effective interactions after integration 
of the phonon and I 2 = 3/4 sector. 

Similarly, exchange interactions in the SO (5) sector are 
also possible: 



H SO{5) = KLt ■ L(0). 



(27) 



Here, K denotes the coupling constants. This can be 
obtained after integrating the sector of I 2 = 0. Note 
that due to the mismatch in the eigenvalue of the spin 
j, a term nj • n(0) cannot exist, when the spin SU(2) or 
time reversal symmetry is present. 

In order to understand that the interaction (|27|) re- 
spects the original symmetry of the Hamiltonian (fl}, 
namely the charge SU(2) symmetry, the spin SU(2) sym- 
metry, and the Z2 symmetry, we expand L/ • L(0) in the 
fermion representations: 



Lj-L(O) = 



i{[/ + ( n (0) St (0)+ P ;(0)p T (0) 

r * [E (4(0)^(0) +pt (ok (o) 

ir 2 {[/ + ( S; (0) St (0)- m (0)p t (0) 

[4(o) SCT (o)- P t ( oK(o)]} 

7 

{[/ + ( w (o> t (o) + . 

[pt(0) Sff (0) + h.c 



h.c. 



)-»]} 



E 



E 



MO)Pt(O) 

} 



h.c. 



h.c. 



+ 



-TV 



]T[4(0K(0)-h.c. 



(28) 



= I • [l(0) + t,I'(0) - T x I"(oj\ - \r v VV{Q). 

(29) 

Here, 1(0) is the local axial charge for the conduction 
electrons. As shown in Appendix [Bl I'(0) is the local 
longitudinal-flavor axial charge, while I"(0) is the local 
transverse-flavor axial charge. All three transform as vec- 
tors under the charge SU(2) rotations, while 2?(0) is in- 
variant; see Appendix [Bl Since I transforms as a vector 
and V is invariant under the charge SU(2) rotations, Eq. 
d2H) is invariant under the charge SU(2) operations. As 
for the parity, I, r 2 , V, 1(0), and I'(0) are even parity, 
while t x , r y , I"(0) and T>(0) are odd parity. Thus, Eq. 
(1291) is even parity, i.e., invariant under the inversion op- 
eration. Finally, since all the terms are spin singlet, Eq. 
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(|29| is invariant under the spin SU(2) operations. These 
facts confirm that L(0) • L/ is invariant under the original 
symmetry. 

In addition to the symmetry of the Hamiltonian (TTJ), 
the exchange interactions (|26|) and (l27l) have an addi- 
tional flavor symmetry. The flavor symmetry operation 
is defined as \ft)++ |JJ-) and s ff f>p„. There are two kinds 
of operators in terms of the flavor symmetry: even or 
odd. Even- flavor operators are denoted by F = 0, while 
odd ones are denoted by F = 1 in Table Q] The values 
of F are related to the eigenvalues of the flavor transfor- 
mation as (— 1) F . As will be investigated in Sec. Ill B 71 
this symmetry breaking drives the system away from the 
critical line. 

Note that there is also flavor SU(2) symmetry in 
Eq. (|27p. when we define the local flavor operators 
as Fj =(— Lf 4 , Lj 4 , Lj 3 ). It is also important to note 
that the phonon flavor operator f cannot directly cou- 
ple with conduction electron flavor F„ in the presence 
of the charge SU(2) symmetry. This means that simple 
flavor-exchange Kondo couplings 

^flavor = JfT ■ F(0) (30) 

never appear under the presence of particle-hole symme- 
try and .ffflavor is absent also away from the critical line. 
In this sense, even away from the critical line, the screen- 
ing processes of impurity degrees of freedom in the model 
(JTJ with the charge SU(2) symmetry are not the same as 
those in the flavor Kondo model such as in the two-level 
systems Away from the critical line with the charge 
SU(2) symmetry, a possible form in the (phonon-only) 
flavor interaction is Ising like: 

^flavor = (J?T X + J)t v )F» <® = (J?T X + J*T„)23(0)(31) 

This is because only D(0) = F v (0) = L 24 (0) is a charge 
and also a spin SU(2) singlet operator in the "local den- 
sity" form constructed by the conduction electron oper- 
ators. Since this is an Ising interaction, when only this 
term is present, nonmagnetic "Kondo effects" never oc- 
cur. 



understand, when we consider the process from large U 
limit, since for large U the relevant operator is expected 
to be the spin S. 

We find that there is an alternative way to derive the 
same NFL spectra, i.e., SO(5) spinor 4 fusion^ *(|JJ. 

),-\1\),-flfl\1\),flfl\^))- Since primary states in the 
SO(5)2 sector are 1, 4, or 5, the fusion rule is 1 — > 4, 
4 — » 1 © 5, and 5 — > 4. These are obtained by discarding 
10 and 16 representations in the SO(5) direct products: 
4<g)4 = l©5©10 and 5 ® 4 = 4 © 16. Indeed, this 
fusion rule generates the same low-energy spectra as the 
spin-1/2 fusion, as shown in Table ITOb) . 

As for double fusions^! of spin 1/2 and SO (5) 4, 
the two different fusions lead to the identical operator 
conten t 34 ! 37 shown in Table ILTl (c). All the operators in 
Table |n](c) that satisfy the SU(2)©SO(5) symmetry can 
be present along the critical line. We expect that, for 
large U, the dominant leading irrelevant operator is in 
the spin sector, while for small U, it is in the SO(5) sec- 
tor, since S is not active for small U. 

Now, one might wonder what the difference between 
the two fusions is. We consider that these two are simply 
equivalent. In order to understand this, we show an ex- 
ample of this kind of situation in an impurity Anderson 
model, when the Coulomb interaction U varies from oo 
to while maintaining particle-hole symmetry. 

As is well known, the ground-state spectra of the 
Anderson model from the strong to the weak coupling 
regime are the same as those in the Kondo model with 
the charge quantum number being shifted. The spectra 
in the Kondo model are obtained by a spin-1/2 fusion via 
direct spin absorption^ Alternatively, the same spectra 
can be obtained by an axial charge i = 1/2 fusion, i.e., 
a 7r/2 phase shift. In the Anderson model, both spin 
and charge degrees of freedom are present. When only 
exchange-type interactions are considered in the context 
of the coarse-grained Hamiltonian, there are two types 
of such interactions: J S S ■ 3(0) and J C I • 1(0). Here J(0) 
and 1(0) represent the spin and axial charge current for 
conduction electrons, respectively. For large U, only the 



5. Fusions 

In this subsection, we introduce two different 
fusion o 30 ' 31 that derive the spectra obtained in the NRG 
calculations along the critical lineJ^r— Although, due to 
the fact that the model ([T]) is not described by a sim- 
ple exchange Hamiltonian, we cannot carry out a direct 
impurity absorption as in exchange models such as multi- 
channel and spin-3/2 multipolar Kondo problem a 30 ! 31 ' 34 , 
we will show that two fusions indeed derive the same NFL 
spectra as in the 2CKM. 

The first candidate of the fusion is spin-1/2 fusion 
in the SU(2)2 sector, which is the same as the case of 
2CKM 3 ^i, leading to the NFL spectra shown in Ta- 
ble [TT] (b). This fusion is physically natural and easy to 



TABLE II: (a) Spectra of the free Hamiltonian for a non- 
degenerate ground state, (b) The spectra at the NFL fixed 
point. The energy Eq and E are measured in the unit of 
ttvf/1. (c) Operator contents at the NFL fixed point. A is 
the scaling dimension of the operators labeled by the quantum 
numbers j and the dimension of the irreducible representation 
in the SO(5) group. 
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sector with I 2 = is relevant. Then, the spin-spin ex- 
change coupling describes the low-energy physics, and 
thus the model reduced to the Kondo model. For small 
U, the charge-charge exchange interaction J c becomes 
compatible with the spin-spin exchange interaction J a . 

Now, one can realize that there are similarities between 
the Anderson model and the present one; J s corresponds 
to J in Eq. (H>D and J c to K and 1(0) (I) to L(0)(L/) 
in Eq. (l27l) . As for the fusion process, the axial charge 
fusion, i.e., the tt/2 phase shift, in the Anderson model 
corresponds to the SO(5)-4 fusion in the present model. 

The spin-1/2 and SO(5)-4 fusions introduced above 
are equivalent in the sense that the spin-1/2 fusion and 
the 7r/2 phase shift are equivalent in the Kondo or An- 
derson model. Our answer to the question "what is go- 
ing on for small U ¥' , which is the main motivation in 
this paper, is that the NFL spectra of 2CKM can be ob- 
tained via the nonmagnetic SO(5)-4 fusion, and thus, we 
can interpret the low-energy physics for small U as the 
nonmagnetic SO (5) Kondo effects in the same way that 
the conventional Kondo effects can be interpreted as the 
strong potential scattering with the phase shift tt/2. Of 
course, one can still interpret it as a magnetic one but 
the SO(5)-4 fusion is much better for understanding the 
physics for small U, since what makes the low-energy 
physics for small U different from that for the Kondo 
regime is the nonmagnetic degrees of freedom. Indeed, 
as will be shown in Sec. IIII1 the residual interactions 
around the NFL fixed point for small U are governed by 
the operators in the SO (5) sector rather than those in 
the spin sector. 

6. Leading irrelevant operators 

Low-temperature thermodynamic properties are gov- 
erned by leading irrelevant operators around the fixed 
point. In the Kondo regime, i.e., for large U, the lead- 
ing irrelevant operator should be in the spin sector, and 
it is J_i • <fi s , with 4> s being spin-SU(2) primary fields 
with the dimension A = 1/2 and the quantum numbers 
(j, SO(5)) = (1, 1) in Table M (c). In total, the dimen- 
sion of this operator is 3/2. In the presence of this oper- 
ator, it is well know n 30 ' 31 that the impurity contribution 
to the specific heat C is proportional to — Tin J 1 , and 
the magnetic susceptibility Xs diverges logarithmically 
% s ~ — InT at low temperatures. 

For small U, we expect that the operator in the spin 
sector does not play an important role, and thus, op- 
erators in the SO(5) sector dominate. Then, the situa- 
tion is analogous to the spin-3/2 dipole-octupole Kondo 
models Since the first descendants of primary fields 4>s 
with (j,SO(5)) = (0,5) and A = 1/2 in Table M (c), 
L_i</>5, cannot form an SO(5) singlet, the leading irrel- 
evant operator is the energy-momentum tensor in the 
SO(5) sector at the impurity site :L(0) • L(0):. The lead- 
ing dimension of this operator is 2, i.e., the "Fermi liq- 
uid" like interaction^ This readily indicates that the 



low-temperature impurity specific heat C cx T. As in- 
vestigated in Ref. 34, the susceptibilities of the SO (5) 
vector 5, <j>&, diverge logarithmically ~ —hxT, indicating 
the divergence of the susceptibility of n/. The suscepti- 
bility of the SO (5) generators is independent of T at low 
temperatures, since the dimension of the (0, 10) operator 
in Table ILT1 (c) is A = 1. Thus, the susceptibilities of Lj 
are Fermi liquid like. We call these behaviors SO(5) NFL 
hereafter. 

Here, we notice that the SO (5) vector in our model 
corresponds to T y I a and TpV with a — x,y, or z and 
[3 = x, or z, see Table U and Eq. (j2"3")l . Physically, these 
operators correspond to bipolaron fluctuations, which are 
coupled fluctuations of the flavor and the axial charge. 
In the original variables b and b\ t x ± ir y are, roughly 
speaking, related to b and bK 

In principle, there exist both terms J_i • <f> s and 
L(0) -L(0) for general values of U, since the original inter- 
action is in complex form of the spin and SO (5) degrees 
of freedom and also terms that cannot be described by 
exchange forms. What varies as a function of param- 
eters, e.g., U along the critical line, is the coefficients 
of these operators in the effective Hamiltonian near the 
fixed point. Such a situation is represented by the fol- 
lowing residual effective Hamiltonian: 

6H e g = A s J_i • <j) s + Az,L_i • L_ 1; (32) 

where A s and Al depend on microscopic parameters such 
as U, Vq, and V\. We have retained the leading term of 
L(0) ■ L(0) in the second term in Eq. (|3"2")l and we have 
not included a term J i • J i, which has the dimension 2, 
since it is sufficient to include only the leading irrelevant 
operators in each of the sectors in the following analysis. 
As we investigated above, the relative magnitude of the 
two terms varies, and the first term is dominant in the 
Kondo regime, while the second one prevails in the weak 
coupling regime. 

An interesting crossover is expected especially in the 
impurity contributions to specific heat C . Since two sec- 
tors are decoupled, C is the sum of the contribution of 
each sector : 29 ' 31 

C = -7.T]n(^) + 1L T. (33) 

Here, T s is a dynamically generated energy scale in the 
spin sector, which is proportional to the Kondo tempera- 
ture in the Kondo regime. The parameter 7 s (7l) is pro- 
portional to A 2 (A_l) i 29 ' 31 As analyzed by Johannesson et 
al.^2. A 2 ~ 1/T S and X L ~ 1/T L , where T L is the "Kondo 
temperature" for the SO(5) sector. The crossover tem- 
perature T* can be defined by the temperature where the 
magnitudes of the two terms in Eq. (|33p are equal, and 
is given by 

T* = T s exp ( - — Y (34) 

For T < T* , the specific heat due to the spin sector dom- 
inates, and thus, C/T diverges logarithmically. However, 
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for sufficiently small U, j l /j s ~ T S /T L > 1, thus, T* is 
never reached in a realistic temperature range and C/T 
stays constant at low temperatures. This is a marked 
contrast between the NFL in the 2CKM and the SO(5) 
NFL. 

Finally, let us comment on the "secondary-diverging" 
susceptibility in each of the parameter regime. Even for 
large U, the susceptibility of ^5, indeed, diverges log- 
arithmically. We call this divergence "secondary", since 
the coefficient of this part is expected to be very small for 
large U; 38 ~ -[hx(T/T L )]/T L with T L ~ U/2. The same 
is true for the spin susceptibility for small U. There, the 
system is in the valence fluctuation regime and the spin 
susceptibility is ~ — [ln(T/T s )]/T s with T s proportional 
to the hybridization width. 



7. Stability of the fixed points 

In this subsection, we investigate the stability of the 
NFL fixed point derived in Sec. Ill B 61 against various 
perturbations. 

First, we investigate symmetry breaking fields. In Ta- 
ble [H] (c), there are SO(5)-5 primary fields 4>5 with the 
dimension A = 1/2. Thus, when the SO (5) symmetry is 
broken, a term h 5 ■ cf> 5 appears in the effective Hamilto- 
nian, which is relevant, and thus the NFL fixed point is 
unstable against this perturbation. Practically speaking, 
the SO (5) symmetry-breaking field /15 includes the inver- 
sion symmetry-breaking field and the flavor (even-odd) 
symmetry-breaking one. 

In the presence of the inversion symmetry breaking 
field, (^5)2 with the quantum number (j,i,P,F) — 
(0,0,1,0) appears in the effective Hamiltonian. Here 
(05 )i represents the ith field in the five component vec- 
tor ^5. With regard to the flavor symmetry breaking 
field, (^5)4 with (j,i,P,F) — (0,0,0,1) appears in the 
effective Hamiltonian. The Hamiltonian ([1]) has inver- 
sion symmetry, while the flavor symmetry is higher than 
that of the original Hamiltonian ([T]) and is realized only 
along the critical line. Thus, (^5)2 cannot appear even 
away from the critical line, while ((^5)4 can appear away 
from the critical line and (^5)4 ((^5)4 is the perturbation 
that makes the NFL fixed points unstable. This flavor 
symmetry breaking causes the energy difference between 
ID") and |J|)- The flavor symmetry breaking is, indeed, 
consistent with the spectra in the NRG and in the small 
clusters as analyzed in Sec. Ill B 31 

Another relevant field is the magnetic field h, since 
there is an SO(5)-singlet and spin-1 primary fields <f> s in 
Table HI (c) with A = 1/2, which couple with h. This is 
the same as in the 2CKM22 and we do not discuss it in 
detail. 

Finally, particle-hole asymmetry is marginal, since the 
operator with (j,i iz ,P,F) — (0, lo,0, 0) is classified in 
[j, SO(5)| = (0,10) and the dimension is A = 1 in Ta- 



ble |TT] (c). This operator breaks the SO(5) symmetry. 
It is well known that the potential scattering VI z (0) is 
absorbed in phase shifts^, and thus, the NFL properties 
are not affected except for the specific heat for small U, as 
we discuss below. When particle-hole symmetry is bro- 
ken by charge-conserved perturbations such as VI z (0), in 
addition to (|32[) with anisotropic exchange interactions 
in the SO(5) sector [see, Eq. ([56")) ] . anisotropic flavor ex- 
change interactions are allowed to appear, leading to the 
additional residual interactions, 

SHtg ~ g x L 3 _\($5)2 + g y L 2 Mfch + 9zL 2 \{$ h )i. (35) 

Here, (— L 3 ^, L 24 -^ L 2 ^) = F_i is the flavor current op- 
erator defined in Sec. Ill B 21 and gi(i = x, y, or z) is con- 
stant proportional to the symmetry breaking field ~ V. 
Note that this form is not SO(5) invariant. However, it 
is still invariant under the inversion and the spin SU(2) 
operations, and also the total charge is conserved. Since 
the scaling dimension of (f3"5|) is A = 3/2, the impurity 
contribution of the specific heat of this term is similar to 
that of the spin sector in Eq. . 

Second, we investigate exchange anisotropies in the 
spin and the SO (5) sectors. The irrelevance of the spin 
exchange anisotropy is explained in the same way as in 
the case of the magnetic 2CKM>2 2 i As for the exchange 
anisotropy in the SO (5) sector, when the anisotropy ex- 
ists, the isotropic effective interaction KL](0) • L(0) is 
replaced, 

ifLj-L(O) -> XL/ • L(0) + SK't z I ■ I'(0) 

-SK"t x I ■ I"(0) - ^SK"'t v VV(0), (36) 

where 5K\ 8K" and 6K'" are the deviations from the 
isotropic interactions. Note that Eq. (|56"|) is charge- 
SU(2) invariant. This anisotropy generates residual in- 
teractions such as L_iL_i with some sets of ab. This 
operator has dimension 2, i.e., it is irrelevant. 

III. NUMERICAL RENORMALIZATION 
GROUP RESULTS 

In this section, we examine the crossover from the 
SO (5) NFL to the NFL in the magnetic 2CKM as U in- 
creases along the critical line by using Wilson's NRG. 23 
One of the advantage in using NRG is that we can ob- 
tain information about the scaling dimensions A of lead- 
ing irrelevant operators around fixed points by analyzing 
spectra obtained in NRG^iS The details of the NRG 
method are explained in a previous paper<21 Here, we 
will analyze variations of NRG spectra along the critical 
line. A detailed analysis of physical quantities will be 
reported elsewhere^ 

In this section, we will show data for three different 
parameters: large U = 0.81), small U = 0.021), and 
intermediate U = 0.21). Here, D is related to half of the 
band width D of conduction electrons for both s and p 
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bands as D — D(A + 1)/(2A), and the Fermi energy is 
at the middle of the band. Here, A is a discretization 
parameter in NRG and we use A = 3. For each value 
of U , we tune the hybridization V\ to realize the NFL 
fixed point, while Vb = 0.2£> and fl = 0.2D are fixed. 
The resulting Vi's are Vi ~ 0.161811? for U = 0.02/3, 
Vl ~ 0.15623.5 for U = 0.21), and V x ~ 0.15473/3 for 
17 = 0.8/3. 

In the calculations, we utilize the spin rotational sym- 
metry to restore the states and 5000 states labeled by the 
set of quantum numbers (j, i z , P) are kept at each itera- 
tion in NRG. As for the number of local phonon states, 
we use 20 phonon states in our calculations. 

A. NRG spectra 




5 10 15 20 25 



N 



Let us show the NRG spectra En as a function of 
the renormalization group step N in Fig. [21 Apart 
from differences in crossover scale iVo(~15, 9, and 3 for 
U/D = 0.02, 0.2, and 0.8, respectively), all three spec- 
tra converge on the NFL spectra 7?nfl shown in Table 
HTl (b). Down to the lowest energy scale, we also confirm 
that the impurity contribution to the entropy is In y/2, as 
expected in the 2CKM 25 , (see, Fig. g}. For the smallest 
U, the crossover step N is large and this is due to the 
fact that for N < No, the nonmagnetic "local moment" 
fixed point is realizedi 21 i 42 

In Table Mil the low-energy states at the NFL 
fixed point and their energy with eigenvalues of the 
z-component of the total axial charge, the total spin, 
and the total parity for U/D = 0.02 are listed. For 
other two i7's, the results are very similar. NRG en- 
ergy eigenvalues and their quantum number are consis- 
tent with those derived by the BCFT. The energy spec- 
tra for the odd- A" NRG step and even-A' one are iden- 
tical within the numerical accuracy except for the par- 
ity eigenvalue P. The eigenvalue P for the even-Af step 
is obtained from that for the odd- A" one simply by in- 
terchanging even (P = 0) and odd (P = 1) for all the 
states. Here, the ground state of the free electron sys- 
tem for the odd (even) A" is non-degenerate (degener- 
ate). This is easily understood by noting that the free 
spectra for even A" can be obtained by shifting the con- 
duction electron charge Q S ( P ) to Q S ( P ) + L 2 ^ where Q s (p) 
is the charge for s(p)-wave conduction electrons. As a 
result, the total parity P = mod(Q p + b^b, 2) is shifted to 
->■ mod(Q p + l + 6t6,2) = mod(P + 1, 2). Thus, the even 
(odd) parity states are simply relabeled as odd (even) 
parity states and then the fusion process leads to the 
NFL spectra with P replaced by mod(P + 1,2). 

The point we address in the following is how En varies 
as a function of A" for the three parameters. Near the 
fixed point, the spectrum at the step N, En, is repre- 
sented as 

En = -Enfl + SEn, (37) 
and the deviation from the fixed point value SEn is given 



FIG. 2: (Color online) NRG spectra as a function of the 
renormalization group step N with N odd for U = 0.813(0), 
U = 0.25(A) and U = 0.02Z)(v). The spectra are uniformly 
scaled such that the first excited energy becomes 1/8 for large 
N, and the identical factor is used for all three. For all three, 
the spectra converge on the NFL spectra /?nfl (dotted lines) 
for large N. The small deviations from -Bnfl are due to trun- 
cation errors in the NRG calculations. 

^ ,,23.39.40 

E(A r -l)JV 
A r A 5 — . (38) 

r 

Here, r identifies leading irrelevant operators appearing 
near the fixed point. For the single channel Kondo model, 
A r = 2. 23 In magnetic 2CKM, A r = 3/2^ which rep- 
resents the "slower" renormalization than in the single 
channel model. Note that A r coincides with the scaling 
dimension of the operator. In the following, we analyze 
SEn in details and examine the crossover predicted in 
Sec. ITTBl 



B. Crossover along the critical line 

The deviation from the fixed point SEn contains in- 
formation about leading irrelevant operators as shown in 
Eq. (f3"8"]) . In order to evaluate SEn, we use En +2 — En'- 

E n+ 2 -E n =J2 KA-^^iA 1 -^ - 1). (39) 

r 

Thus, when there is a dominant leading irrelevant 

term with A r = A, En+2 — En is proportional to 
A -(A-i)JV/2 _ SEn Here; we uge ( N + 2). and iv-step 

eigenvalues in Eq. (|39[) . since, in NRG, there is even-odd 
alternation in the spectra. 

Figure [3] shows \E N+2 — E N \ of the low-energy states 
for the three parameters. Each state is labeled by the 
spin and the SO(5) indices: (j, SO(5)). For the largest 
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TABLE III: Spectra E N at the NFL fixed point for odd TV = 
29 (even N = 30) and U = 0.0273. The parity P and the 
NRG spectra En for even TV are shown in parentheses. Em 
is scaled such that the first excited energy with * becomes 
1/8 in order to compare En with the spectra obtained by the 
BCFT 75bcft. 



P SO (5) Sbcft E 29 (E 3 



0(1 



0(1 
1(0 
0(1 
1(0 



-1 





1 

_ 1 

J 

2 



1(0 
0(1 
1(0 
1(0 
1(0 



0(1 
1(0 
0(1 
1(0 



0(1 



-1 
-1 
-1 






1 
1 
1 



0(1 
0(1 
1(0 
0(1 
0(1 
1(0 
1(0 
0(1 
0(1 
1(0 



0(1 



(0) 



0.12503 
0.12500 
0.12503 
0.125 



(0.125*) 
(0.12504) 
(0.12500) 
(0.12504) 



0.50155 
0.50154 
0.50124 
0.50155 
0.50155 



(0.50154) 
(0.50153) 
(0.50136) 
(0.50154) 
(0.50154) 



0.63098 
0.63100 
0.63098 
0.63100 



(0.63104) 
(0.63100) 
(0.63104) 
(0.63100) 



1 



1.00306 (1.00304) 
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1 1.02581 
1.02616 
1.02615 
1.02581 
1.02616 
1.02615 
1.02580 
1.02581 
1.02616 
1.02615 



(1.02600) 
(1.02620) 
(1.02619) 
(1.02600) 
(1.02620) 
(1.02619) 
(1.02599) 
(1.02600) 
(1.02620) 
(1.02619) 



1 



1.03984 (1.03992) 



U = 0.873, it is clear that the scaling dimension A is 
A = 3/2, and thus, the NFL is described by the magnetic 
2CKM as investigated in Sec. Ill Bl For the smallest U — 
0.021?, the dimension of leading irrelevant operator is 
A = 2, since for most of the states the N dependence 
is ~ A _JV / 2 . This is consistent with our analysis in Sec. 
Ill B 61 The Ar N / 2 dependence is due to the existence of 
the nonmagnetic SO (5) residual interaction L_i -L_i. In 
principle, there is a possibility that the term J_i • J_i 
is the origin of the A _JV / 2 dependence. This possibility, 
however, is unlikely from a physical standpoint. It is 
unphysical that only the coefficient of the leading term 
in the spin sector is suppressed, while that of the sub- 
leading term in the same sector is not, as U decreases. 

One may find that some of the curves follow the A~ N ^ 4 
dependence for U — 0.0273, but the absolute value is 
very small, i.e., |A|| is very small, where we use A r as 
the index r. Although, in principle, there exist contribu- 
tions of matrix elements of the operators in the effective 
Hamiltonian to |7?/v+2 — it is unlikely that the small 
absolute value is only due to the matrix elements, and 
thus, we neglect them in the following analysis. This N 
dependence, A _Ar / 4 , must originate in the magnetic in- 
teractions J_i • 4> s , since in the SO(5) sector there are no 
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FIG. 3: (Color online) |75jv+2 — 7Jjv| for 15 low-energy states 
as a function of the renormalization group step N with N 
odd for (a) U = 0.873, (b) U = 0.273, and (c) U = 0.0273. 
Each state is labeled by the spin j and the dimension of the 
irreducible representation in the SO(5) group. A straight line 
~ A~ N ^ 2 is characteristic of the irrelevant operator in the 
SO(5) sector L-i -L— i, while ~ A _JV//4 is characteristic of the 



operator in the spin sector J_ 
300 low-energy states. 



Inset: the same plot for 



such operators that generate this A^ dependence. Note 
that in NRG the step Af is related to the energy scale 
as DA~ N / 2 , and, for example, A^ = 20 corresponds to 



3 -iO£, ^ i 6 g x iq-5 D At TV = 20, the absolute value 
of IA2I is more than 100 times larger than | A3 |. Thus, we 
expect that there the spin degrees of freedom have no ef- 
fect on, for example, the specific heat since | A2 1 >• |A||. 
As for the spin susceptibility, a logarithmic divergence 
with a very small coefficient is expected, reflecting the 
small A 3 . This is similar to the case of the flavor suscep- 
tibility in the two-channel Anderson models 

As expected, the situation changes as U increases. For 
U = 0.2D, it is clear that both A2 and A3 are present 
with similar magnitudes. Around N ~ 22, the crossover 
from the SO (5) NFL to the NFL in 2CKM occurs. Thus, 
from our NRG calculations, it is clear that the profile of 
the leading irrelevant operators changes smoothly from 
the weak-coupling regime to the Kondo regime. These 
results confirm the results in Sec. IIIBI 

Finally, we discuss the impurity contributions of spe- 
cific heat C and the impurity entropy S. Figure 0] shows 
the temperature dependence of {S — Sq)/T and S for 
the three parameters of U with So ~ In \/2- Since 
(S — Sq)/T ~ C/T at low temperatures, it represents 
C/T for T < Tq. Crossover temperatures are defined as 
T = Dh~ N °/ 2 . For large U/D = 0.8, C/T at low tem- 
peratures diverges logarithmically and this is consistent 
with the conventional magnetic 2CKM. As expected from 
the results of the scaling dimension of leading irrelevant 
operators, the temperature dependence of C/T changes 
with decreasing U. One can clearly see that it is constant 
at low temperatures for U/D = 0.02 and 0.2. For T < T , 
the temperature dependence of entropy for U/D = 0.02 
and 0.2 is well described by a single-scale function of 
(T/T ) ~ (T/T L ), while that for U/D = 0.8 has a differ- 
ent functional form, since To ~ T s for U/D = 0.8. For 
U/D = 0.2, the temperature dependence changes from 
~ const, to — InT, and the crossover temperature has 
been defined as T* in Eq. (|34|). This is the crossover 
from the SO(5)-operator dominant regime to the spin- 
operator dominant regime. Below T ~ T*, the tempera- 
ture dependence of S is not described by the single-scale 
function of (T/Tl), although the deviation in S is very 
small ~ — TlnT. Note that the impurity entropy S is 
S ~ ln-s/2 for the temperature where C/T = const, ap- 
pears. These results confirm our BCFT predictions and 
the importance of nonmagnetic SO(5) degrees of freedom 
for small U. 



IV. DISCUSSION AND SUMMARY 

In this paper, we have analyzed low-energy critical 
theory in a two-channel Anderson model with phonon- 
assisted hybridization on the basis of BCFT and NRG. 
One important finding is that nonmagnetic SO (5) degrees 
of freedom are constructed in the I <g> t sector, which are 
"hidden" in the 2CKM due to the Hilbert space restric- 
tion, and also the conduction electron part is rewritten 
by the SO(5) currents. These nonmagnetic degrees of 
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FIG. 4: (Color online) Temperature dependence of (S — So)/T 
for the three values of U. Temperature T is scaled by Tq; 
Tq/D = 0.192 for U/D = 0.8, T /D = 0.00713 for U/D = 0.2, 
and Tq/D = 0.000264 for U/D = 0.02. Inset: Temperature 
dependence of S. 



freedom are important for small U. We have demon- 
strated that SO(5)-4 fusion gives exactly the same NFL 
spectra as those in the magnetic 2CKM. The difference 
between the spin-1/2 and the SO(5)-4 fusions have been 
discussed and we interpret the fusion as simply equiva- 
lent ones, noticing that the spin-1/2 fusion and the 7r/2 
phase shift in the Anderson model arc equivalent. A full 
understanding of the exact fusion process will require a 
more sophisticated analysis. 

In the form of residual effective interaction (|32p . the 
crossover between small and large U can be described by 
changes in the coefficients X s and Xl in Eq. (|32p . Note 
that the residual interactions in the SO (5) sector never 
appear if we map the model to the magnetic 2CKM, since 
there is no A^L(0) • L(0) term in the perturbation ex- 
pansion in the spin sector. Using physical intuition, we 
predict that the SO (5) sector is more important than the 
spin sector for small U, leading to linear specific heat at 
low temperatures. This has been checked by the NRG 
calculations; the scaling dimension of leading irrelevant 
operators varies from A = 3/2toA = 2as[/ decreases, 
and the impurity contribution of specific heat is propor- 
tional to temperature for small U. 

The difference between the present NFL for small U 
and the NFL in the two-level systems should be noted, 
although they both have a nonmagnetic origin. It is 
well known that the NFL spectra in the two-level Kondo 
model is derived via flavor-1/2 fusion^ The spectra is 
the same as those in the magnetic 2CKM when the spin 
and the flavor sector are interchanged, while the present 
NFL spectra are exactly the same as those in the mag- 
netic 2CKM. The scaling dimension of the leading irrele- 
vant operator in the two-level Kondo model is 3/2, which 
leads to a logarithmic diverging specific heat coefficient. 
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Thus, nonmagnetic Kondo phenomena in the two models 
should be distinguished and the microscopic mechanism 
for the NFL in each of the models is quite different, i.e., 
flavor SU(2) and nonmagnetic SO(5) exchange interac- 
tions. 

Our BCFT analysis can also be applicable to the an- 
harmonic model investigated in Ref. 24, since even in 
the anharmonic phonon model, parity P is a good quan- 
tum number and the phonon states in the effective the- 
ory would be described by r as in a similar manner to 
that of the present analysis. With regard to the gen- 
eralization of this model to the more realistic one, it is 
interesting to take into account optical modes that couple 
with localized electrons. This electron-phonon coupling 
reduces the bare Coulomb repulsion. When it is suffi- 
ciently large, the effective Coulomb interaction becomes 
attractive, and thus, it is possible to realize another NFL 
fixed point in which the spin and charge sector are inter- 
changed from the present NFL for U > 0. With regard 
to a lattice generalization of the present model, it is also 
interesting to examine whether some composite pairing 
operators listed in Table U condensate into exotic super- 
conducting states. 

In summary, we have investigated the microscopic ori- 
gin of the line of non- Fermi liquid fixed points found pre- 
viously by numerical simulations We have succeeded 
in constructing nonmagnetic SO (5) degrees of freedom, 
and, on the basis of boundary conformal field theory, we 
have pointed out that, for the weak electron-electron cor- 
relation regime, the non-Fermi liquid can be interpreted 
as an SO(5) non-Fermi liquid, which crosses overs to a 
non-Fermi liquid in the Kondo regime. We have also an- 
alyzed the difference in the leading irrelevant operators 
as U varies, and indeed we have confirmed it by numeri- 
cal simulations. In particular, the impurity contributions 
to the specific heat are proportional to temperature T 
for small U, while they are proportional to — TlnT for 
large U. The present results demonstrate that it is im- 
portant to take into account not only single degrees of 
freedom, e.g., only a phonon, but also complex degrees 
of freedom formed both by electrons and phonons in the 
Kondo problems in electron-phonon coupled systems. 



Appendix A: Matrix representations of SO(5) 
Generators and Vectors 



In this Appendix, we summarize the definitions of 
SO(5) matrices. We follow the convention used by Wu, 
et alM 



1. Generators: 10 representation 



The SO(5) generators C ab define all the SO(5) rotation 
and are given by L ab = T ab /2 with 



-i 12 



V2:i 



V25 



'ooob 

10 
10 

a o o 0/ 

/ —1 ON 
1 
-10 

V o i o 0/ 

/l 
0-100 
10 

\0 -1, 



il3 



/ 



-i 

-i 



/° 




Vi 



^^ 

1 

0j 






-1 







-1 

0/ 



p24 _ 



r 34 = 



^45 





-i 

i 

v-i 0/ 

(-\ 0\ 
0-100 
10 

V o o o i/ 

i 0\ 

-i 

i 

-i 0/ 

/0 1 
10 
1 

Vo o i 



0/ 







I 
-i 



1 



(Al) 



£ ab, s are ten-dimensional adjoint representation 10 and 
satisfy the following SO (5) commutation relations: 
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[C ab ,C cd ] = -i(5 bc £ ad - S ac C bd - 6 bd C ac + 6 ad C bc ), (A2) 
= if ab ' cd ^f £ ef , (A3) 



where the repeated indices are assumed to be summed 
over and l<a<6<5, l<c<d<5, and 1 < e < / < 
5. C ab with a > b should be regarded as C ab = —£ ba on 
the right-hand side of Eq. (|A"2j) . 
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2. Vectors: 5 representations 

The 5 representation is an S0(5) vector and is repre- 
sented by the following five matrices: 



r 1 



r d = 



r 4 = 



o o i o\ 
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-i 0/ 

'0 -i 0\ 
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v -i Oj 

/0 1 0\ 
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10 
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In terms of r a 's, r af "s are represented as 
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(A4) 



[r a ,r 6 ]. 



(A5) 



The commutation relations between T a and r fcc are given 
by 

[r\T bc ] = -2i(sgii(c~a)S ab T c + S gii(b-a)6 ac T b ).(A6) 



Appendix B: Axial Charge 

In this Appendix, we summarize various spherical ten- 
sors with respect to the axial charge. They appear as a 
part of the SO(5) degrees of freedom in the main text. 



1. Spherical tensors of axial charge symmetry 

Spherical tensor operators Tm in the axial charge sec- 
tor arc defined by 



(Bl) 



[W.T®] - V(lTm)(l±m + l)T«> ±1 . (B2) 

Here, I to t is defined by Eqs. ([3]) and (|4]) and I is an 
integer and is called the rank of operator Tm and \m\ < I. 
Operators with I = are scalar, i.e., invariant under 
the charge SU(2) operations, while operators with I > 1 
transform as rank-^ tensors. In particular, operators with 
/ = 1 transform as vectors. A trivial example is the axial 
charge of conduction electrons I(x), 

I+(x) = s\{x)s\(x)+p\(x)p\{x) = [I-{x)]\ (B4) 



I(x) denotes the contributions of the conduction 
electrons to Itotr^ and the rank-1 tensor with 
[TlV , T (1) , T^] = [/_ (x)/V2, h (x), -1+ (x)/V2] . Thus, 
the quantum number of I(a;) is (Jj z ,i,P) = (0o,l,0), 
where, is the eigenvalue of spin (axial charge) with 
the z-component j z and P is the parity. Here, the quan- 
tum numbers in the spin sector can be determined in the 
same way as in the axial charge sector. The quantum 
numbers of the /-electron axial charge I are the same as 
those of I(a;). In the following, we will list various spher- 
ical tensors that appear in the SO (5) degrees of freedom 
analyzed in the main text. 



2. Longitudinal-flavor axial charge 

We define longitudinal-flavor axial charge I'(x), which 
is related to the SO(5) generator L 12 (x), L 23 (x) and 



L (x), and appears in Eq. (|2"5j) as 

40*0 = \j2[ s °(xW(x)-pUx)Pa(x)], (B5) 

l' + (x) = s\(x)s\(x) - p\(x)p\(x) = (x)]l (B6) 

The parity of these operators is even, since p a appears 
as quadratic forms in Eqs. (|B5I) and (|B6|) . In terms 
of the spherical tensors, l'(x) is the rank 1 tensor with 
[T^It£\t[ X) ]=[I'_{x)/^/2, I' z {x), -P + {x)/V2]. Thus, 
the quantum number of I'(x) is (jj zl i,P) — (0q, 1,0). 



3. Transverse-flavor axial charge 

We define transverse-flavor axial charge l"{x), which 
is related to L 14 (x), L 34 (x) and L A5 (x), and appears in 
Eq. U2SJ as 

W*) = \ E [pU^W(x) + S \(x)p a {x)] , (B7) 

(T 

I'l{x) = s\(x)p\(x)+p\(x)sl(x) = [l'l(x)]K (B8) 

The parity of these operators is odd, since I"(x) 
includes one p a in each term. In terms of the 
spherical tensors, l"(x) is the rank 1 tensor and 
[TiV , T (1) , Tf»] = [I» (x)/V2, 1'J(x) , (x)/V2] . Thus, 
the quantum number of I" (a;) is (jj z ,i,P) — (0q, 1,1). 



4. Flavor singlet 

We define a flavor singlet operator T>(x) = L 24 (x), 
which appears in Eq. (I2TR) as 

a 

The parity of this operator is odd, as is evident from 
the fact that T>{x) includes one p a in each term. Since 
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[I to t± : V(x)] = and [I totz ,'D(x)) = 0, T>(x) is a scalar 
operator and the quantum number of T>(x) is ,i,P) = 
(00,0,1). 



5. Longitudinal spin-flavor singlet 

We define a longitudinal spin-flavor singlet operator 
T>'(x) = n A (x), i.e., the fourth component of the SO(5) 
vector n(x) in Table Q] as 

= ±ya[sUx) S(T (x)-pUx)p a (x)\.(B10) 

CT 

The parity of this operator is even, as is evident from the 
fact that T>'{x) includes zero or two p CT 's in each term. 
Since [I tot± ,V'(x)\ = and [I totz ,V (x)} = 0, V'(x) is 
a scalar operator and the quantum number of T>'(x) is 
(j^,i,P) = (lo,0,0). 



6. Transverse spin-flavor singlet 

We define a transverse spin-flavor singlet operator 
V"(x) = n 2 (x) in Table Has 

WW = \y^[pl{x)s a (x) + sl{x) Pa (x)\.{mi) 

The parity of this operator is odd, since T>"(x) includes 
one p a in each term. Since [itot±, T^"( x )\ = and 
[It t Z iV" {x)\ — 0, V"{x) is a scalar operator and the 
quantum number of V" (x) is (jj z ,i,P) — (lo, 0, 1). 



7. Transverse-spin-flavor axial charge 

We define transverse-spin-flavor axial charge I"'(x), 
which is related to n (x), n 3 (x) and n 5 (x) in Table U 
as 

iy, 

2 ^ 



(B12) 



TO = l -y<T[ P l{x)s a {x)-sl{x)p„{x) 

(7 

I™{x) = -i[s\(x)p\(x) + sl(x)p\(x)] =[/ w (.x)]t.(B13) 

The parity of these operators is odd, since I"'(x) 
includes one p a in each term. In terms of the 
spherical tensors, V"(x) is the rank 1 tensor with 

[TlV, TV]=[I>!!(x)/V2, I'»(x), -I'i'(x)/V2}. Thus, 
the quantum number of l"'(x) is (jj z ,i,P) = (1q, 1, 1). 



Appendix C: Primary states of SO(5)2 Kac-Moody 
algebra 

In this Appendix, we briefly show that the primary 
states for k = 2 SO(5) Kac-Moody algebra are 1, 4, and 
5 representations. 

Since the rank of SO (5) group is 2 and thus the Cartan 
subalgebra consists of Hi = Lj 5 and H 2 = £q 3 , we can 



label primary states by their eigenvalues, hi and h 2 , and 
denote them as \hi, /12) with hi and h 2 being integers or 
half-integers. Since , ^2) is primary, L^lhi,^) = 
for n > 0. Ladder operators are defined as 



(i) 



J 



Ji 2 > = i 



iL 



24 



[J. 



Wit 



(CI) 

= [jf]t,(C2) 



and they satisfy 

[Hi,f 1) }=0. 
[Hi,J^} = - 



[H 2 ,jW] = -J 



(i) 



J 



(2) 



(2)1 



J 



(2) 



(C3) 
(C4) 



The commutation relations Eqs. (|C3I) and (IC4[) indicate 
J_ /12) 1^1, h 2 — 1) and ji 2 ^ \h±, h 2 ) oz\hi — 1, /12 + 



Similarly, we can define another set of ladder operators 



as 



J 



(i) - 



Li\ +i Li\=[m, 



(C5) 



J 



(2) _ 



L 3 + \ - Z# - - iL» ^[j| 2) ]t(C6) 



where L"* = [LljJ^. Straightforward calculations show 



that they satisfy 
and also satisfy 



(i) 



-J. 



(2) 



[^2,ji 2) 



2ifo 



fW 



(2) 



and 



[J- 



Hi-H 2 - k/2. 



(C7) 
(C8) 

(C9) 
(C10) 



Now, let us consider the norm of descendant states. Since 
the norm is positive, we obtain 



Ijj 1 ^!,^)! 2 = (/i 1 ,/i 2 |Ji 1J J| 1J |/ ll ,/i 2 )= 

= {hiM[j { -\j+ ] ]\hi,h 2 ), 

= -2h 2 + k >0, (Gil) 

where at the second line we have used L a t 1 \h-i,h< 2 ) = 
and at the third line, Eq. (|C9|) and (hi, h 2 \hi, h 2 ) = 1 

~(2) 

have been used. Similar calculations for J\_ \hi,h 2 ) lead 
to 



W 7W1 



h 2 -hi + k/2 > 0. 



(C12) 



It is clear that irreducible representations with the large 
dimension cannot satisfy Eqs. (|C11|) and (IC12I) . since, in 
general, such irreducible representations have large \hi\ 
and \h 2 \. Indeed, Eqs. (|C11I) and (|C12|) indicate that 
there are three irreducible representations, and they are 
the identity 1, the spinor 4, and the vector 5. Thus, 
primary states in the SO(5)2 sector belong to 1, 4, or 5 
representations . 
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